The angular spin current and its physical consequences 
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5-H . We find that in order to completely describe the spin transport, apart from spin current (or linear 

spin current), one has to introduce the angular spin current. The two spin currents respectively 
describe the translational and rotational motion of a spin. The definitions of these spin current 
densities are given and their physical properties are discussed. Both spin current densities appear 
naturally in the spin continuity equation. Moreover we predict that the angular spin current can 
also induce an electric field E, and in particular E scales as 1/r 2 at large distance r, whereas the E 
field generated from the linear spin current goes as 1/r 3 . 



• PACS numbers: 72.25.-b, 85.75.-d, 73.23.-b 

Recently, a new sub-discipline of condensed matter physics, spintronics, is emerging rapidly and generating great 
interests. |lj, |2( The spin current, the most important physical quantity in spintronics, has been extensively studied. 
Many interesting and fundamental phenomena, such as the spin Hall effect [^, 0, Q and the spin precession 
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systems with spin-orbit coupling, have been discovered and are under further study. 

As for the charge current, the definition of the local charge current density j e (r,t) = Re[w(r, t)v9(r, £)] and its 
q ■ continuity equation Jjp(r, t) + V • j e (r, t) = is well-known in physics. Here ^(r, t) is the electronic wave function, 
v = r is the velocity operator, and p(r,t) = ^tvp j s the charge density. This continuity equation is the consequence 
' of charge invariance, i.e. when an electron moves from one place to another, its charge remains the same. However, 
' in the spin transport, there are still a lot of debates over what is the correct definition for spin current. The problem 
7-H ' stems from the lack of understanding of the spin continuity equation. Recently, some studies have begun investigation 
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, in this direction, e.g. a semi-classical description of the spin continuity equation has been proposed [81 13 ■ 
, In this paper, we study the quantum version of the continuity equation and the definition of local spin current 
density. We find that in order to completely describe the spin transport, one has to use two quantities, the linear spin 
current and the angular spin current, which respectively describe the translational and rotational motion of a spin. 
The conventional linear spin current has been extensively studied. However, the angular spin current has not been 
investigated before. The definition of two spin current densities are given and they appear naturally in the quantum 
spin continuity equation. Moreover, we predict that the angular spin current can generate an electric field similar as 
with the linear spin current, and thus contains physical consequences. 
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Before studying the spin current in a quantum system, we first consider the classical case. Consider a classical 
particle having a vector m (e.g. the classical magnetic moment, etc.) with its magnitude |m| fixed under the particle 
motion. To completely describe this vector flow (see Fig. la), we need three quantities: the local vector density 
M(r, t) = p(r, t)m(r, t), the velocity v(r, t), and the angular velocity (3(r, t). Here p(v, t) is the particle density, and v 
and to describe the translational and rotational motion, respectively. Since |m| is a constant, the continuity equation 
is (see Fig. la): fffy MdV = — § s dS • vM + ffj v ujx MdV, or it can be rewritten in the differential form: 

-M(r, t) = -V • v(r, t)M(r, t) + tD(r, t) x M(r, <), (1) 

where vM is a tensor, and its element (vM)ij = mMj. 
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Now we study the electronic spin s in the quantum case. Consider an arbitrary wave function ^(r, t). The local 
spin density s at the position r and time t is: s(r, t) = ^(r, t)s^(r, t), where s — |<r with 3 being the Pauli matrices. 
The time-derivative of s(r, t) is: 



d K 
—s(r,t) = - 
dt v ' ' 2 



dt 



+ ^ r 3^\ . (2) 



dt 

From the Schrodinger equation, we have ^^(r, t) — -^rH^i(r,t) and ^^(r, i) = -^j-[H 1 ii(r,t)] J< . Notice here the 
transposition in the symbol f only acts on the spin indexes. By using the above two equations, the Eq.(2) changes 
into: 

(d/dt)s(r, t) = [^3H^ - (H^3^}/2i. (3) 
In the derivation below, we use the following Hamiltonian: 

H=^-+V(r)+3*B+^z*(3xp). (4) 
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Note that our results are independent of this specific choice of the Hamiltonian. [lOj In Eq. (4) the 1st and 2nd terms 
are the kinetic energy and potential energy. The 3rd term is the Zeeman energy due to a magnetic field, and the 
last term is the Rashba spin-orbit coupling , which has been extensively studied recently H 0, Q 

Next we 

substitute the Hamiltonian Eq.(4) into Eq.(3), and calculate the right side of Eq.(3) one term by one term: 

(1) . For the 1st term of Eq.(4), |^: 

&dr(jf /2m)V - \(f/2m)^3^ 

= -(h 2 /2m){&ffV • V* - [V • W]^*} 
= -(K 2 /2m)V • {^V*?* - [V*]^*} 

= -iKV • Re{^(p/m)3$>} (5) 

Here the dotted multiply between the vector and the tensor is defined as usually, e.g. (A • B)j = J2i AiBij with the 
vector A and the tensor B. 

(2) . For the 2nd term of Eq.(4), V^r): its corresponding value is zero. 

(3) . For the 3rd term, 3 • B: 

¥3(3 • B)V - [3 • m\ r 3^ 

= ¥ (33 • B - B»33)^ = 2i^B x <?* (6) 

(4) . For the last term, j-z • (3 x p) = %3 • (px z): 

^33 • (px z)^ -[3 • (px z)^3^ 
= ¥(p x z) • 33$ + ^2i(p x z) x a* — [3 • (px £)\E , ] t a\l/ 
= -iK{¥(V • (z x 3)3^ + [V*]* »(zx 3)3^} + 2i¥(p x z) x 3^ 
= -iKV • {4^(1 x 3)3^} + 2i¥(p x z) x <r* 

In the similar method, one also has: 

^33 • (px z)^ - [3 • (p x z)^3^ 
= -iKV • {¥(z x 3)3^}* + 2i[(p x z) x 3^]^ 
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Therefore, we obtain the result of the last term: 



(a/h){^aa i(px z)^ - [a • (px z)^^} 

2ia 



iaV • Re^Hz x ct)<t<I/ H Re^Up x z) x ct* (7) 

ft 



To summarize the above four terms, the Eq.(3) changes into: 



MftfB+^x I] x (8) 
ft 



Introducing a tensor j s (r, and a vector j'^r, f): 



j.(r,t) = fle{*t[Z + "(£x50]l^} (9) 
to ft 

J w (r,t) = i2e{*tH[J+ ^xf)]xj*}, (10) 
ft ft 



then Eq.(3) reduces to: 

d 



^s(r,t)=-V.j 5 (r,i)+i a) (r,t), (11) 



or it can also be rewritten in the integral form: 



d_ 



sdV = - <b dS.j s + j^dV. (12) 
v Js JJJv 



Due to the fact that v — = £ + f(z x a) and ^<r = ^-[ct, if] = + § p x z] x <r, Eqs.(9) and (10) become: 

j a (r,t) = fle{*t(r,t)tffo(r,t)} (13) 
jL(r,f) = i?e{^ t (d?/dt)*} = fle-f* 1 ^ x <r$<}, (14) 



where Q = f (px 

Eq.(ll) is the spin continuity equation, which is very similar with the classic vector continuity equation (1). In fact, 
this spin continuity equation (11) is the consequence of invariance of the spin magnitude \s\, i.e. when an electron 
makes a motion, either translation or rotation, its spin magnitude \s\ = \ remains a constant. The two quantities 
j s (r, t) and j^ir^t) in Eq.(ll), which are defined in Eqs.(13,14) respectively, describe the translational and rotational 
motion of a spin at the location r and the time t. They will be named the linear and the angular spin current densities 
accordingly 1 13J. In fact, the linear spin current j s (r,t) is identical with the conventional spin current investigated in 
recent studies|l4j. However, the angular spin current J w (r, f) is new and has not been investigated before. |13j 

Next, we discuss certain properties of j s (r, t) and j w (r, t). Notice that j w (r, t) which describe the rotational motion 
of the spin plays a parallel role in comparison with the conventional linear spin current j s (r, t) for the spin transport. 
(1) The linear spin current is a tensor. Its element, e.g. j s ,xy, represents an electron moving along the x direction 
with its spin in the y direction (see Fig. lb). The angular spin current j u is a vector. In Fig.lc, its element j u x 
describes the rotational motion of the spin in the y direction and the angular velocity lu in the — z direction. (2) Why 
is it necessary to introduce two quantities instead of one to describe a spin motion? The reason is the same as for 
completely describing the motion of a classic vector or a rigid body. There two quantities are used: the linear velocity 
for the translational motion and the angular velocity for the rotation. It is similar for the present quantum case, both 
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j s (r, t) and ju>(j, t) are required to describe the motion of a quantum spin. (3) From the linear spin current density 
j s (r, t), one can calculate (or say how much) the linear spin current I s flowing through a surface S (see Fig.le): 
if = JJ S dS • j s . However, the behavior for the angular spin current is different. From the density j u (r,t), it is 
meaningless to determine how much the angular spin current flowing through a surface S, because the angular spin 
current describes the rotational motion not the movement. On the other hand, one can calculate the total angular 
spin current T% in a volume V from j u ; T% — Jff v j^(r,t)dV. (4) If the system is in a steady state, j s and j u are 
independent of the time t, and 4zs(r,t) = 0. Then the spin continuity equation (11) reduces to: V • j s = j u or 
§ s dS • j s = JJJy jujdV . This means that the total linear spin current flowing out of a closed surface is equal to the 
total angular spin current enclosed. If to further consider a quasi one dimensional (ID) system (see Fig.le), then 
one has if — if = ijf. (5) It is easy to prove that the spin currents in the present definitions of Eqs.(13,14) are 
invariant under a space coordinate transformation. (6) The linear spin current density j s = Re{ty^ vsty} gives both 
the spin direction and the direction of spin movement, so it completely describes the translational motion. However, 
the angular spin current density, j u = -Rej^' 4| 1 ^ , } = Re{ty^uj x sty} involves the vector product of uj x s, not the 
tensor ujs. Is it correct or sufficient to describe the rotational motion? For example, the rotational motion of Fig.le 
with the spin s in the y direction and the angular velocity uj in the — z direction is different from the one in Fig. Id in 
which s is in the z direction and uj is in the y direction, but their angular spin currents are completely the same. Shall 
we distinguish them? It turns out that the physical results produced by the above two rotational motions (Fig.le 
and Id) are indeed the same. For instance, the induced electric field by them is identical since a spin s has only the 
direction but no size (see detail discussion below). Thus, the vector j u is sufficient to describe the rotational motion, 
and no tensor is necessary. 

Recently, theoretic studies have suggested that the (linear) spin current can induce an electric field Can 
the angular spin current also induce an electric field? If so this gives a way of detecting the angular spin current. 
Following, we study this question by using the method of equivalent magnetic chargejl^- Let us consider a steady- 
state angular spin current element jujdV at the origin. Associated with the spin s, there is a magnetic moment (MM) 
fh = g[LsS = 2g ^ B s where \ib is the Bohr magneton. Thus, corresponding to j w , there is also a angular MM current 
jmujdV = 2<3 f L B jujdV . From above discussions, we already know that j mul (or j u ) comes from the rotational motion of a 
MM m (or s) (see Fig.le and Id), and j mu = Q x fh (or j u — uj x s). Under the method of equivalent magnetic charge, 
the MM m is equivalent to two magnetic charges: one with magnetic charge +q located at 8n m and the other with — q 
at —6n rn (see Fig. If). h m is the unit vector of fh and S is a tiny length. The angular MM current j muJ is equivalent 
to two magnetic charge currents: one is j+ q — hjq5\uj\ sin# at the location Sh m , the other is j_ g = njg<5|a;| sin# at 
—5h m (see Fig. If), with hj being the unit vector of j mu and 8 the angle between uj and fh. In our previous work,^^ 
we have given the formulae of the electric field induced by a magnetic charge current. The electric field induced by 
jmujdV can be calculated by adding the contributions from the two magnetic charge currents. Let 5 — > 0, and note 
that 2qS — ► |m| and |w||m| sin0 = |j mw |, we obtain the electric field E u generated by an element of the angular spin 
current j^dV: 

E u = — — / ; = / judV x — (15) 



4:71 J r 3 h 

We also rewrite the electric field E s generated by an element of the linear spin current using the tensor j s :Q| 

E s = Vx isdVm-^, (16) 



Below we emphasize three points: (i) In the large r case, the electric field E u decays as 1/r 2 . Note that the field from 
a linear spin current E s goes as 1/r 3 . In fact, in terms of generating an electric field, the angular spin current is as 
effective as a magnetic charge current, (ii) In the steady-state case, the total electric field Et = E u + E s contains 
the property: § c Et * dl = 0, where C is an arbitrary close contour not passing through the region of spin current. 
However, for each E^ or E s , § c E^, • dl or § C E S • dl can be non-zero, (iii) As mentioned above, a angular spin 
current j w may consist of different Co and s (see Fig.lc and Id). However, the resulting electric field only depends 
on j u — Co x s. This is because a spin vector contains only a direction and a magnitude, but not a spatial size (i.e. 
the distance <5 approaches to zero). In the limit 5 — > 0, both magnetic charge currents j± q reduce to Co x to/2 at the 
origin. Therefore, the overall effect of the rotational motion is only related to to x to, not separately on Co and to. 
Hence it is enough to describe the spin rotational motion by using a vector uj X s, instead of a tensor Cos. 

Now we give an example of applying the above formulas, Eqs.(13,14), to calculate the spin currents and the electric 
fields generated by them. Let us consider a quasi ID quantum wire having the Rashba spin orbit coupling, and its 
Hamiltonian is: 

H = ^ + V(y,z) + ~ a zPx , (17) 



Here the other Rashba term — j^<J x p z is neglected because z-direction is quantized jg. Let $ be a stationary wave 
function 

*(*) = — (yd**,)* j <Piv,z), (18) 

where kji = am/h 2 and <p(y, z) is the bound state wave function in the confined y and z directions, ^(r) represents 
the spinprecession in the x-y plane while moving along the x axis and this wave function was discussed in recent 
studies. |y, 13 Using Eqs.(13,14), the spin current densities of the wave function ^(r) are easily obtained. There are 
only two non-zero elements of j s (r): 

iaxx(r) = (h 2 k/2m)\y(y,z)\ 2 cos2k R x 1 (19) 
j sxy (r) = {h 2 k/2m)\y(y,z)\ 2 sm2k R x. (20) 

The non-zero elements of j w (r) are: 

3um(T) = -{h 2 kk R /m)\^(y, z)\ 2 sm2k R x, (21) 
j uy (r) = (h 2 kk R /m)\y{y,z)\ 2 cos 2k R x. (22) 

Those spin current densities confirm with the intuitive picture of an electron motion, precession in the x-y plane and 
movement in the x direction (see Fig. 2a). Moreover, they indeed satisfy —V • j s (r) + j w (r) = 0. 

Let us calculate the induced electric fields at the location r = (x,y, z) by the above spin currents. Substituting the 
spin currents of Eqs. (19-22) into Eqs. (15,16) and assuming the transverse sizes of the ID wire are much smaller than 



y 2 + z 2 , the induced fields E w and E s can be obtained as 

p v in\ ■ ~ ( zcosx I y K {(3) 

Euj — —aKi(p)smx 1 



E- - ay sini (o, ^ + zF((3), - yF(, 1 1 
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where i = 2fcjji (i = x,y, z), (3 = \/ y 1 + z 1 , F(f3) — 1 + ^ , the constant a — 2po91 l B Kkk\p s /~Km, p s is the linear 
density of moving electrons under the bias of an external voltage, and Kq and K\ are the Bessel functions. The total 
electric field is: Et = E u + E s = —^^\7{ zsl ^ x Ki(f3)}. The total electric field Et represents the one generated by a 
ID wire of electric dipole moment p e = (0, 0, c since) at the x axis (see Fig. 2b), where c is a constant. It is obvious that 
V x Et = 0, i.e. § c Et • dl = 0. However, in general § c E u • dl and § c E s • dl are separately non-zero. Finally, we 
estimate the magnitude of Et- We use parameters consistent with realistic experimental samples. Take the Rashba 
parameter a = 3 x 10 _11 eVm (corresponding to fc^=l/100nm for m — 0.036m e ), p s — 10 6 /m (i.e. one moving 
electron per 1000 nm in length), and k = hp = 10 8 /m. The electric potential difference between the two points A and 
B (see Fig. 2b) is about O.OlpV, where the positions of A and B are 2^(1,0,0.01) and ^(^f, 0, 0.01). This value 
of the potential is measurable with today's technology. Furthermore, with the above parameters the electric field Et 
at A or B is about lOV/m which is rather large. 

In summary, we have introduced the angular spin current and its role in the spin continuity equation. We point 
out that the angular spin current can also induce an electric field and its E field scales as 1/r 2 at large r. 
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FIG. 1: (Color online) (a) Schematic diagram for a classic vector flow, (b) The linear spin current element j s ,xy (c) and (d) 
The angular spin current element j u ,x- (e) The spin current in a quasi ID quantum wire, (f) The currents of two magnetic 
charges that are equivalent to a angular MM current. 

FIG. 2: (Color online) (a) Schematic diagram for the spin precession in the x-y plane while moving along the x axis, (b) A 
ID wire of electric dipole moment p e . This configuration will generate an electric field equivalent to the field from the spin 
currents in (a). 
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